Effect of next-nearest neighbor hopping on the spin dynamics in antiferromagnets 
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Recently, inelastic neutron scattering (INS) experiments on the insulating parent compounds of 
high-r c materials were analyzed to extract the value of the superexchange constant J. Starting 
point of the analysis was the nearest-neighbor Heisenberg model. Motivated by recent ARPES 
experiments, we consider the effects of a next-nearest neighbor hopping, t' in the strong coupling 
limit of the spin-density wave formalism, where it leads to an antiferromagnetic exchange J' > 
between next-nearest neighbor spins. We show that a finite J' (a) changes the spin-wave velocity and 
thus the value of J extracted from the INS data, yielding almost identical values for J in La2Cu04 
and YBa2Cu30e.i5, and (b) leads to a considerable change in the high-frequency spin-dynamics. 
Finally, we calculate the 1/S quantum corrections in the presence of J' and discuss their relevance 
for the INS data. 
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The strong evidence for antiferromagnetic fluctuations 
as the source of the pairing mechanism in the high- 
T c materials [Qj2) makes it important to understand 
the variation of the energy scale J of these fluctua- 
tions in different compounds. The materials for which 
J can be easily determined are the insulating parent 
compounds of the high-T c materials. Here, the energy 
scale J of the spin-wave mode can be extracted from 
Raman JMJ or inelastic neutron scattering (INS) experi- 
ments (jJT Recently, INS experiments on the insu- 
lating parent compounds La2Cu04, Nd2Cu04, Pr2Cu04 
and YBa2Cu306.i5, which possess an antiferromagnetic 
ground state, were analyzed using the Heisenberg model 
with nearest-neighbor exchange J. Combining the exper- 
imentally measured spin-wave velocity c sw with the the- 
oretical expression c sw — \[2J of the Heisenberg model, 
values of J were obtained which range from J rs 125 meV 
in YBa 2 Cu30 6 .i 5 to J w 160 meV in Nd 2 Cu0 4 . How- 
ever, various theoretical studies ^|J^] of angle-resolved 
photoemission (ARPES) experiments || as well as band- 
structure calculations B suggest, that the experimen- 
tally measured fcrmionic dispersion provides evidence for 
a finite next-nearest neighbor hopping, t', with t'/t f» 
-0.25 for La 2 Cu0 4 and t'/t m -0.45 for YBa 2 Cu 3 6 .i5. 
In the strong coupling limit of the insulating parent com- 
pounds, t 1 leads to an antiferromagnetic coupling J 1 > 
between next-nearest neighbors spins which changes c sw 
and in turn the value of J extracted from the INS data. 

In this communication we investigate the effect of J' on 
the spin dynamics in single and double-layer compounds. 
We show that J' not only changes the spin-wave velocity 
but also the spin-dynamics at higher frequencies. We re- 
analyze the INS data and demonstrate that J extracted 
in the presence of a finite J' is actually larger than previ- 
ously assumed Q. In particular, we obtain almost identi- 
cal values for J in La2Cu04 and YBa2Cu306.is. Further- 
more, we compute 1/S (quantum) corrections [jlO 11 for 
various physical quantities. These corrections increase 
with increasing J', which is expected since J' > frus- 



trates the system and should thus lead to an increase of 
quantum fluctuations. Finally, we discuss our results in 
the context of the experimentally observed strong renor- 
malization of the transverse spin susceptibility. 

We study the effects of t' on the spin dynamics by 
computing the transverse susceptibility within the spin- 
density wave (SDW) formalism at T = 0. Since this 
formalism is thoroughly described in the literature Jl2| , 
we will outline it only briefly. For simplicity we first con- 
sider the case of a single-layer compound. The starting 
point for the description of the Cu02 planes in the in- 
sulating parent compounds is the effective 2D one-band 
Hubbard model |||l3| a ^ na lf filling, given by 



(1) 



Here a is the spin index and t% j is the hopping inte- 
gral which we assume to act between nearest and next- 
nearest neighbors (t and t', respectively). In the strong 
coupling limit, t <C U the ground state of the Hub- 
bard model is antifcrromagnctically ordered, and one can 
use the antiferromagnetic long range order parameter 
(S z ) = X)fe( c l c fe+Q) t° decouple the interaction term in 
Eq. ([l]) . Diagonalizing then the quadratic form by means 
of a unitary transformation one obtains two electronic 
bands for the conduction and valence fermions whose 
mean-field energy dispersions are given by 



Et = ± 



\/W) 2 + A 2 + 



k ' 



where ef = (e k ± e k+Q )/2, A = U(S Z ), and 

gfe = — 2t(cos k x + cos k y ) — At 1 cos k x cos k y — fj, 



(2) 



(3) 



Here E? ,v is the dispersion of the conduction and va- 
lence fermions, respectively, e k is the dispersion of free 
fermions, and [i is the chemical potential. 

Next we calculate the transverse susceptibility which 
is defined by 



1 



X + -(q,q',t)^(TS+(t)SZ q ,(0)) 

where the spin operators are given by 
1 



S„ = 



q - 2 2s C k+q,a CT ":/3 c k,/3 



(4) 



(5) 



and a are the Pauli matrices. Within the RPA approx- 
imation is given by an infinite series of fermionic 
bubble diagrams 12[ , with fermions obeying the disper- 
sion given in Eq.(2^7 Since the susceptibilities with zero 
momentum transfer (i.e. q' = q) and with momentum 
transfer Q = (tt, 7t) (i.e. q' = q + Q) are finite, the eval- 
uation of x -1 reduces to a (2 x 2) algebraic equation. 
One obtains after some tedious algebra (for S=l/2) 



-2J 



x (q,q»w) = — 
x + ~(q,q + Q,w) 



l-Vq- 



- cos(q x )cos(q y ) 



uj 2 — uj 2 (q) + iS 



uj 1 — uj 2 (q) + iS' 



(6) 



where J = 4i /U, J' = 4(t') /U are the superexchange 
constants and Vq — (cos(q x ) + cos(q y ))/2. The magnon 
dispersion w(q) is given by 



w(q) = 2J 



J' 

1 - — (1 - ^(g^)^^)) 



(7) 



A comparison of Eq.(j?]) with the mean- field magnon dis- 
persion of the Heisenberg model with next-nearest neigh- 
bor exchange J' shows, as expected, that in the strong 
coupling limit t 1 leads to an antiferromagnetic coupling 
J' > 0. 

From Eq.(|) we see that the terms ~ J' vanish in the 
limit q = 0, which means that the static, uniform sus- 
ceptibility x± — 1/8J is not affected by J'. In contrast, 
J' changes the spin- wave velocity c sw , which is given by 



(8) 



Using the hydrodynamical relation p s = c sw x±_ we obtain 
for the spin stiffness p s = J(l — 2J'/J)/4 which was 
previously reported in Ref . |jj . 

We now discuss our results in the light of the avail- 
able INS data. In Fig. [I] we present a fit of the magnon 
dispersion Eq.(j^) for t' = — 0.3i (solid line) to the exper- 
imental data on La2Cu04 from Ref. H]. This fit yields 
J = 166 (±6) meV, a value approximately 6% larger than 
in Ref. g, and J' = 15 ± 1 meV. For comparison, we 
also present a fit of the magnon dispersion with t' = 
(dashed line). Note, that only the functional form of the 
high-frequency part of the dispersion is changed by J' . 
Here, J 1 removes the degeneracy along the boundary of 
the magnetic Brillouin zone (MBZ) ((0, 7f) to (it, 0) direc- 
tion), and gives rise to a local maximum along the (0, tt) 
to (jr, 7r) direction (indicated by an arrow). It is, how- 
ever, questionable whether INS experiments are sensitive 
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FIG. 1. Fit to the experimental magnon dispersion for 
t' = -0.3i (solid line) and t' = -0.0 (dashed line). The 
experimental data are taken from Ref. pM. 



enough to detect this difference since the effect is only of 
the order of 10%, and thus might very well lie within the 
experimental errorbars. 

Next, we analyze the local, or q-integrated imaginary 
part of the susceptibility, which is defined by 



^2 J d 2 qlm$ ] x + (q,q,w)j . 



In the low- frequency limit, we obtain 



X» = 



y/2c sw yJl-2J'/J 



l + O 



(9) 



(10) 



Note that on the mean-field level, J' increases the overall 
scale of x"( w ) (this is opposite to the effect of quantum 
corrections, as we will discuss later). In Fig. |^ we present 
a fit of x"(u) with t' = —0.3 to the experimental data 
from Ref. For this comparison we multiplied Eqs.(|^) 
and (10) with Zdg 2 p.%, where g is the Lande- factor, hb is 
the Bohr magneton and Zd represents the overall renor- 
malization due to quantum fluctuations [^,^1). The in- 
set shows the systematic evolution of x"^) with increas- 
ing t'/t. For t' = 0, x"(w) is logarithmically divergent 
at the maximum magnon energy Lo max due to the de- 
generacy of the dispersion along the boundary of the 
MBZ. A finite t' removes this degeneracy (see Fig. [l]) and 
transforms the divergence into a maximum which shifts 
to lower energies with increasing t'/t. We can identify 
three characteristic features of x"(o;), which in addition 
to c SWl provide enough information to extract J' from 
the experimental data. These features are the position 
of the maximum, the cut-off energy, and the position of 
the jump, denoted by an arrow in Fig. | The cut-off en- 
ergy is obviously given by the maximum magnon energy 
^max = 2 J — 2 J', but for small J' yields the same in- 
formation as c sw . To determine the position of the jump 
and the peak we observe that at higher frequencies x"(w) 
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FIG. 2. Fit of the local susceptibility with t' = -0.3t to 
the experimental data from Ref. [^J. The inset shows the 
systematic evolution of x"( w ) with increasing t'/t 



reflects the density of states (DOS). It is then easy to 
check that the jump in x"( w ) arises from the excitation 
of magnons at q = (0, tt) with energy w q = 2 J — 4 J'. 
Similarly, the peak in x"^) is determined by the flat 
maximum of the dispersion along the (0,7r) to (7t, 7t) di- 
rection with energy 



^peak 



2J 



J-2J' 
J + 2J' 



(11) 



However, a fit of x"( w ) to the experimental data (Fig. ||) 
shows that the accessible energies are still about 75meV 
below any of these pronounced features, which makes it 
difficult to extract a precise value for J'. 

We now turn to the analysis of the two-layer systems. 
The only change in our starting point Eq.(Q) is the in- 
troduction of a hopping term t± between the planes. Di- 
agonalizing the Hamiltonian as before, one obtains two 
valence and two conduction bands, whose dispersion is 
given by Eq. (g) with ej^ replaced by ±t±. The sub- 
sequent calculation of the transverse susceptibilities is 
performed along the same lines as before, and has been 
reported earlier for t' = |L5[ ]. Due to the coupling of 
the planes, one obtains a finite out-of-plane susceptibility 
with the two spins in Eq.(|) on different layers. The final 
result for the in-plane and out-of-plane susceptibilities in 
RPA approximation is given by 



xt a (q, q,w) 



j 



1 - J/ q - - c °s(qx)cos(q y )) 



J 



T- 



1-u, 



1 T 2J .7 



(1 - cos(q x )cos(q y )) 



- ^ac(q) + i6 



x?; ct (q,q + Q,w) 



i 



w 2 -^o pt (q)+^ 



± 



l 

^ 2 - w a C (q) + iS 



(12) 



where a — 1 , 2 (upper and lower sign) , denotes the in- 
plane and out-of-plane susceptibilities, respectively. Note 
that x(q, q+ Q, w) changes sign under the exchange of 
layers. The optical and acoustic magnon modes are dou- 
bly degenerate with dispersion 



w ac ,o P t(q) = 2J< 



J±_ J' 



1/2 



1 ± z/ q - —(1 - cos{q x )cos{q v )) 



(13) 



where the upper (lower) sign corresponds to the acoustic 
(optical) mode. The gap A = 2^/JJ± of the optical (or 
even) mode at Q = (jr, tt) is not changed by J', however, 
the spin-wave velocity of the acoustic (or odd) mode is 
changed and now given by 



(14) 



It is important to note that since both A and c sw are ex- 
perimentally measured, J' does not only lead to a renor- 
malization of J, but also of J±. Using Eq.(|l4|) one ob- 
tains (in the limit of small Jj_) 



J± 



V2 



1 - 2- 



A\ 2 
~2 



(15) 



Using t' = — 0.45i, we obtain from a fit of Eq.(|T3|) to 
the experimental data Q, J = 163 meV, J' = 33 meV, 
and Jj_ = 8.4 meV. The value for J is about 30% larger, 
and the value for J± about 25% smaller than the ones 
reported in Ref. It is interesting that the different 
values for t' in La2Cu04 and YBa2Cu30g.i5 lead to al- 
most identical values for the nearest-neighbor exchange 
constants, in contrast to the results in Ref. 0]. 

In Fig. H we present a fit of the local susceptibility in 
the even and odd channel, x'even^dd^) = Xu ± Xi2> to 
the experimental data. The fit in both channels is quite 
good, though the scattering of the data is rather large. A 
finite t' again transforms the divergence of x"( w ) a t u max 
into a maximum. Note, that the two maxima in Xeven odd 
do not occur at the same frequen cy, but are shi fted by a 
finite amount Acu = 2J±(J'/J)/y/l - (2J'/J) 2 w 3meV. 
As before in the case of a single-layer, one can use the 
position of the jump and the maximum together with the 
value for c sw to extract J'. However, as a comparison of 
x" for J' = (dashed line) shows, the experimental data 
have to be extended to higher frequencies before one can 
extract J'. 

Furthermore, using the data in Figs. ^| and we 
compared the local in-plane susceptibility Xii if") i n 
YBa2Cu306.i5 and La2Cu04. We obtain that for uj < 
A = 72 meV, Xii( w ) m YBa 2 Cu306.i5 is only half as 
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FIG. 3. The local susceptibility as a function of uj for the 
even channel and (see inset) the odd channel. The solid line 
corresponds to t'/i = —0.45, whereas the dashed line shows 
the result for t' /t — 0. The experimental data are taken from 
Ref. fll. 



large as in La2Cu04, which is expected since the optical 
mode is gaped out. For 72 meV < w, the optical magnons 
in YBa2Cu30g.i5 can be excited and x'd m both mate- 
rials is approximately the same. 

Finally, we consider quantum corrections in the pres- 
ence of J'. A convenient way to calculate these correc- 
tions is to perform a systematic l/S'-expansion in the 
Heisenberg model, as was done for the case J' = in 
Rcf. Jllj . Considering a finite J' we performed calcula- 
tions analogous to jLlJ. To lowest order in 1/S, i.e., on 
the mean- field level, we recover as expected Eqs.(^|) and 
(0). The next order 1/5 corrections, which appear due 
to the normal ordering of quartic operator terms yield 
the following renormalizations (for J' / J = 0.09) 



2S 

. , 0.188N 



i 



0.123 



J', 



9 S 



Comparing these results with the case J' = 0, we find 
that the spin-wave velocity is only slightly increased 
whereas Zd is decreased by about 10%. Can these re- 
sults explain the strong renormalization of \" observed 
experimentally ? A fit of x"(0) = Zdg 2 ^ 2 3 /V2c sw to the 
experimental data Q showed, that the overall renormal- 
ization Z d — 0.36 of the dynamical spin susceptibility 
in La2Cu04 was much larger than the theoretical value 
Z d = 0.607 obtained in the case J' = @. For J' ^ 
the overall renormalization of x" is given by a classical 



contribution, 1/yl — 2J'/J, and the quantum correc- 
tions Zd, which means that Z e , xp should be compared 



measurable. With J' /J — 0.09, which neglecting small 
corrections to order (1/S) 2 JllJ implies J'/ J = 0.126, we 
obtain Z' = 0.57. Though this value is only about 6% 
smaller than the one in Ref. jnj , it shows that a finite J' 
leads to an improved agreement with the experimental 
data. 

In conclusion we have provided a detailed study of the 
effects of t', and consequently J', on the spin dynamics 
in single and double-layer systems. We have shown that 
J' changes the value for J and J± extracted from the 
experimental data, yielding two almost identical values 
for J in YBa2Cu306.i5 and La2Cu04. We demonstrated 
that the changes due to J' in the high frequency struc- 
ture of provide the possibility to extract values 
for the exchange constants. Precise results, however, will 
require a higher accuracy of the INS data and their ex- 
tension to higher frequencies. Furthermore, we obtained 
that J' leads to stronger quantum corrections. Finally, 
it was recently pointed out that the 17 NMR signal in 
Ni doped YBa2Cu306+rr is sensitive to the momentum 
structure of x'Cl) @h an effect which might be used to 
detect changes in x'{o) due to J'. Further quantitative 
studies in this direction are clearly called for. 
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with Z' = Zd/ y 1 — 2 J'/ J. Here we use the renormal- 
ized J, J' and c sw , since only these are experimentally 



[1] D. Pines, Z. Phys. B 103, 129 (1997). 

[2] D. J. Scalapino, Phys. Rep. 250, 329 (1995). 

[3] K.B. Lyons, P.A. Fleury, L.T. Schneemeyer, and J.V. 

Waszczak, Phys. Rev. Lett. 6 0, 732 (1988). 
[4] S. M. Hayden et al, preprint, cond-mat 9710181. 
[5] P. Bourges, H. Casalta, A.S. Ivanov. and D. Petitgrand, 

Phys. Rev. Lett. 79, 4906 (1997). 
[6] D. Duffy and A. Moreo, Phys. Rev. B 52, 15607 (1995). 
[7] A.V. Chubukov and K. Musaelian, J. Phys. Cond. Mat. 

7, 133 (1995). 
[8] B.O. Wells et al, Phys. Rev. Lett. 74, 964 (1995). 
[9] M. S. Hybertsen, E. B. Stechel, M. Schluter, and D. R. 

Jennison, Phys. Rev. B 41, 11068 (1990). 
[10] R.R.P. Singh, Phys. Rev. B 39, R9760 (1989). 
[11] J. Igarashi, Phys. Rev. B 46, 10763 (1992). 
[12] J.R.Schrieffer, X.G.Wen, and S.C.Zhang, Phys. Rev. B 

39, 11663 (1989); A.V.Chubukov and D.M.Frenkel, Phys. 

Rev. B 46, 11884 (1992). 
[13] PW. Anderson and J.R. Schrieffer, Phys. Today 44, 55 

(1991). 

[14] S. M. Hayden et al, Phys. Rev. Lett. 67, 3622 (1991). 

[15] D.K. Morr, A.V. Chubukov, A. Kampf, and G. Blum- 
berg, Phys. Rev. B 54, 3468 (1996). 

[16] J. Bobroff et al, Phys.Rev. Lett. 79, 2117 (1997); D. 
Morr et al, to appear in Phys. Rev. Lett. 



4 



